A density-functional study of charge doping in WO3 
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The addition of electron donors to the vacant A site of defect-perovskite structure tungsten trioxide 
causes a series of structural and chemical phase transitions; for instance, in the well-known case of 
the sodium tungsten bronzes (Na^WOs). Here we calculate the effect of the addition of electronic 
charge to WO3 without the complication of also including sodium ions. Our density functional 
theory method enables isolation of electronic effects from the additional size, chemical, and disorder 
effects present in experimental samples. Our calculated low-temperature phase diagram between 
1=0 and a;=l moves from the initial low-temperature monoclinic phase through a second (centered) 
monoclinic phase, an orthorhombic phase, a tetragonal antiferroelectric, and an aristotypic cubic 
phase, in broad agreement with the experimentally-observed transformations in Na^WOs. Our work 
confirms that the observed structural transformations are driven primarily by electronic factors. We 
find that the dominant electronic effect is the covalent interaction between the tungsten 5d and 
oxygen 2p orbitals. 

PACS numbers: 61.72. Ww, 71.15.Mb, 71.20.Nr, 71.30.Bh, 71.38.Ht 
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I. INTRODUCTION 

Tungsten trioxide, WO3, is a technologically- 
significant ceramic: it is electrochromioi, and the pos- 
sibility of ion intercalation and dcintcrcalation gives rise 
to several potential applications in devices (for instance, 
as a cathode of rechargable batteries-). In particu- 
lar, it is often used in optical applications due to the 
fact that the color can be changed by doping with elec- 
trons. (This was first observed in 1815 by Berzelius^ in 
H^WOs). Bulk stoichiometric WO3 is yellow-green in 
hue, but the sodium-doped tungsten bronzes, Na^WOs, 
exhibit most colors of the visible spectrum on varying Na 
concentration 4 . In addition, the optical absorbance and 
reflectivity of the material can be modulated by injec- 
tion or extraction of electrons and ions, giving excellent 
control over and tunability to the optical properties^. 

In addition to strongly affecting the color of the ma- 
terial, the incorporation of electron donating ions also 
has a strong effect on the structure. For example, com- 
plex structural behaviour is observed on doping with elec- 
tron donors, such as Na, Li or H. Such dopants occupy 
the vacant perovskite A site and therefore they can be 
easily introduced over a wide range of concentrations, 
from trace concentrations in the so-called e phase up to 
the limiting case of the tungsten bronze, NaWC>3. In 
NaWC>3, there is one donor per formula unit, and an aris- 
totypic perovskite structure is adopted^. Its structure is 
the same as that of Re03, with both solids retaining the 
ideal cubic perovskite structure at all temperatures. In- 
deed, assuming that in NaWC>3 the valence electron of 
Na is donated into the 5rf-type conduction/antibonding 
band of WO3, the two materials are isoelectronic. In 
contrast, WO3 displays both off-centering of the W ion 
from its ideal centrosymmetric position, and Glazeri- 
type tilting transitions. Studies on intermediate concen- 
trations of Na^WO^ suggest a complex phenomenology, 
involving a consistent reduction in the degree of polyhe- 



dral tilting with increasing dopant concentration. No- 
tably, Clarke observed that the room-temperature struc- 
tures of Na K W0 3 (0.62 < x < 0.94) are slightly dis- 
torted from the high-symmetry Pm3m aristotype, and 
proposed a preliminary phase diagram to explain the ob- 
served diffraction data™. 

There have been several prior computational studies 
on WO3, carried out at various levels of theory. Here we 
summarize the results of the density functional theory 
(DFT) studies from the literature. Perhaps the most rel- 
evant for our study is the work by Cora et aliSiiS on the 
electronic structure of cubic WO3, Re03 and NaW03. 
Using the full potential linear muffin tin orbital method, 
with the local density approximation, they showed that 
the band structures of the three compounds are very sim- 
ilar, with the extra electron in ReC>3 and NaWC>3 occupy- 
ing the antibonding conduction band and decreasing the 
metal-oxygen bonding strength. They also showed that 
displacement of the transition metal from its centrosym- 
metric position towards an oxygen ion causes increased 
hybridization of the transition metal 5d t2 g orbitals at 
the bottom of the conduction band with the O 2p orb- 
tials at the top of the valence band. This in turn lowers 
the energy of the top of the valence band, and raises that 
of the bottom of the conduction band, explaining why 
metallic Re03 and NaW03 remain cubic, whereas "d°" 
WOgii has an off-centering distortion. Hjelm et &\H 
used the same method to establish that, in L1WO3 and 
NaW03, rigid band filling of the WO3 conduction band 
occurs, whereas in HWO3 the hydrogens form hydroxide 
units with the oxygen atoms and change the electronic 
structure, de Wijs et ali£ calculated the electronic prop- 
erties of the various experimental structural phases of 
WO3 using a plane wave ultra-soft pscudopotential im- 
plementation of DFT within both the generalized gradi- 
ent and local density approximations. They found that 
increases in pressure are readily accomodated by tilting 
of the octahedra, explaining the small experimental bulk 
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TABLE I: Cutoff radii, r(Q, for the basis sets corresponding 
to W and O species; as the basis set is double-^ polarized, 
there are two basis functions (and therefore radii) per hypo- 
thetical atomic orbital, as specified by quantum numbers n 
and I. All radii are given in bohr. The starred entries are 
polarization orbitals. 

modulus. Also, they observed large W displacements, ac- 
companied by strong rehybridization and changes in the 
electronic band gaps. This latter observation is consis- 
tent with the anomalously large Born effective charges 
calculated for cubic WO3 by Detraux et aliii 

The work presented here builds on these earlier theo- 
retical studies by treating fractional doping effects (i.e., 
the effect of electron concentrations between and 1 per 
W ion) in a fully self-consistent manner. Our compu- 
tational approach allows the isolation of the effects of 
additional valence electrons from other factors, such as 
the presence of donor cations and structural disorder at 
intermediate concentrations. Indeed, it has already been 
shown that simple model calculations including only elec- 
tronic effects can reproduce some aspects of the observed 
structural behavior—. Our main result is that electronic 
effects are able to account fully for the experimentally 
observed structural phase transitions. 



II. METHODOLOGY 

Calculations were performed using the Siesta 
implementatioiii£ii&, of density-functional theory 17 
(DFT), within the Perdew-Zunger 18 parametrization 
of the local-density approximation 1 ^. Core electrons 
were replaced by norm-conserving Trouiller-MartinsSS, 
pseudopotentials, factorised as prescribed by Kleinman 
and BylanderiSi The valence electrons were taken to be 
2s 2 2p 4 for O, and 6s 2 5d 4 for W. Wave-functions were 
expanded in a basis of numerical atomic orbitals of finite 
support 22 . A double-C polarised (DZP) basis set was 
used, whereby two basis functions per valence atomic 
orbital are included, plus extra shells of higher angular 
momentum to allow polarization. This size of basis has 
been found, in prior studies, to give results of similar 
quality to typical plane-wave basis sets in other major 
codes. Details for the pseudopotentials and the basis 
sets are given in Tables III and ITT1 

For fc-point sampling, a cutoff of 10 A was used24; this 
gave 24 independent fc-points in the first Brillouin zone 
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TABLE II: Confinement radii for each of the s, p, d, and / 
channels of our W and O pseudopotentials. All radii are given 
in bohr. The W pseudopotential used also included a partial- 
core correction, with a radius of 1.30 bohr, according to the 
scheme of Louie et al-i 



of the low-temperature monoclinic phase, corresponding 
to a 4x3x3 mesh. This is equivalent to 6 independent 
fc-points within the first Brillouin zone of a 64-atom su- 
percell of the same structure, corresponding to a 2x2x3 
mesh (if one were to neglect degeneracy). Both fc-point 
meshes were generated by the method of Monkhorst and 
Pack 2 ^. The fineness of the real-space grid used for nu- 
meric integration was set to correspond to an energy 
cutoff of 200 Ry. A grid-cell sampling 1 - 6 of four points 
(arranged on a face-centered cubic lattice) was used to 
reduce the space inhomogeneity introduced by the finite 
grid. All calculations were performed using variable-cell 
relaxation, with the convergence criteria set to corre- 
spond to a maximum residual stress of 0.01 GPa, and 
maximum residual force component of 0.04 eV/A. 

In order to test the quality of the pseudopotential and 
basis set used, we first performed calculations on the 
same phases as investigated by de Wijs et al. Results 
of comparable quality to theirs were obtained (see ta- 
ble llllf) . although our results tend to give the expected 
underestimation of the experimental lattice parameters 
whereas theirs, even in the LDA, tend to overestimate. 

In order to investigate the effect of additional dopant 
charge on the bonding character and structure of WO3, 
a 2 x 2 x 1 supercell of the low-temperature monoclinic (e) 
phase was used (the aforementioned 64-atom supercell). 
This contains sixteen W centers, and can be considered to 
be a 2\/2 x 2-^/2 x 2 supercell of the aristotypic perovskite 
structure. As such, it is sufficiently large to encompass 
all ground-state symmetries experimentally observed in 
pure WO3. Additional charge was added by the injection 
of extra electrons into the system, with charge neutrality 
over space being ensured by a corresponding homoge- 
neous positively-charged background. 



III. CALCULATED PHASE DIAGRAM 

We performed calculations for numbers of dopant elec- 
trons between 1 and 16 per supercell, (corresponding to 
electron concentrations in the range < x < 1 per W 
center) in single-electron steps. Full, variable-cell, struc- 
tural relaxations were performed for each doping level. A 
series of five phases was observed over the doping range. 
Here we describe the five structures in turn, starting with 
the highest doping level. We refer to figure ^ which 
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TABLE III: Lattice parameters, experimental, from our sim- 
ulations, and from the work of de Wija-^ for (hypothetical) 
cubic, low temperature monoclinic, room temperature mon- 
oclinic, high-temperature tetragonal and room temperature 
triclinic phases of WO3. All distances are given in A. Note 
that de Wijs' calculations were performed under the local- 
density approximation for monoclinic phases, and under the 
generalized gradient approximation for triclinic and tetrago- 
nal phases. 
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shows (a) our calculated lattice parameters (c' is the c 
lattice parameter multiplied by y2 in order to normalize 
it with respect to the a and b parameters) , (b) our calcu- 
lated unit cell volumes, and (c) our W-0 bond lengths. 
Our suggested phase boundaries are also shown in fig- 
ure H(a) as solid vertical lines. 

At maximal doping, we obtain the perovskite aristo- 
typic (i.e. the highest symmetry phase in the series: in 
the case of perovskites, primitive cubic, a = b = c, one 
formula unit per unit cell) structure, with space group 
Pm3m. The three lattice parameters arc equal, as are 
the W-0 bond lengths. 

With decreasing dopant concentration the unit cell 
volume decreases slightly, until at a: « yi, we find 
a symmetry-breaking transition to an antiferroelectric 
tetragonal phase with PA/nmm symmetry. The W atom 
moves off the center of its Oq octahedron in the [001] 
direction [as shown in figure |5] (a)], resulting in differ- 
ent W-0 bond lengths and a lowering of the space group 
of the system. Interestingly, all three lattice parameters 
continue to have the same effective length with respect to 
our supercell (within the accuracy of our computations) 
in spite of the inequivalence of one of them due to the 
reduction in symmetry. This can be ascribed to the fact 
that the oxygen framework distorts very little in these 
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FIG. 1: (a) Lattice parameters (for normalized cubic supercell 
with respect to the e phase) of WO3; (b) Volumes of 2x2x2 
supercell of WO3; (c) Bond lengths within WO6 octahedra, 
all versus electron dopant concentration. 



transitions, and it is the oxygen framework that is largely 
responsible for the volume and shape of the unit cell. 
The antiferroelectric tetragonal phase persists down to a 
dopant concentration of approximately i electron per W, 
with a gradual increase in off-centering with decreasing 
dopant concentration. By x — h, the off-centering of the 
W along the [001] direction has increased to around 0.08 

A. 

The antiferroelectric tetragonal phase becomes unsta- 
ble under i electrons per W, and an orthorhombic phase 
is stabilised. There is also a clear discontinuity in the 
volume per unit cell at this doping concentration. This 
arises from a tilt of the WOq octahedra towards [110]- 
type directions which in turn causes the loss of the 
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FIG. 2: Offcentering in W0 6 octahedra; (a) along [001], (b) 
along both [001] and [110]-type directions 



four-fold rotation in the c direction, lowering the space 
group to Pnma. Initially, the c' and a lattice param- 
eters (within our pseudocubic supercell), despite being 
symmetrically inequivalent, remain essentially equal in 
length; this can be ascribed to the absence of distortion 
of the WC>6 octahedra, which originally tilt (around an 
axis parallel to b) as a rigid unit without substantial de- 
formation. These parameters become unequal at x w X ; 
however, we do not find a lowering of the space group 
caused by this. 

Below x w |, a monoclinic P2\/c phase (similar 
to that identified by de Wijsi^), related to the low- 
temperature e phase in WO3, is stable. The phase bound- 
ary from the orthorhombic to monoclinic phase is char- 
acterized by an increase in the W-0 displacement along 
[001] (seen in the increased bond splitting in figure^M) 
as well as a marked rotation around z (loosely speaking in 
the x—y plane). The Glazer tilt is expressible as a~b~c~ . 
The monoclinic symmetry persists through the remain- 
der of the phase diagram; however an additional phase 
boundary exists at x i. Here an additional splitting 
in the W-0 bond lengths indicates offcentering of the W 
atom in a [110]-type direction (as shown in figure[3 (b)) 
and there is also a second discontinuity in the cell vol- 
ume. This corresponds to a loss of screw axes along b, 
and hence a further lowering of the space group to Pc, 
the e phase. 

Order parameters can be defined for the transitions 
between each of these phases as shown in figure [3J The 
order parameter for the Pc to P2\/c transition (figure |3| 
(a) )is the difference between the W-0 bondlengths in 
the x — y plane (xy\ and xy2)- The order parameter for 
the P2\/c to Pnma transition is, as expected, the devia- 
tion of the (3 angle from 90 ; this goes sharply to zero at 
x « I (figure El(h)). The order parameter for the Pnma 
to PA/nmm phase is the length difference between the 
two remaining unequal lattice constants, \b — a\ (figure |3| 
(c)). Finally, the ideal cubic Pto3to perovskite structure 
is reached when the W moves to its centrosymmetric po- 
sition, and the order parameter for the tetragonal to the 
cubic phase is the magnitude of the W off-center displace- 
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FIG. 3: Order parameters for the doping-induced phase tran- 
sitions in WO3. (a) Pc - P2i/c, in A, where xyi and xyi are 
the inequivalent W-0 bond lengths in the the x — y plane; (b) 
P2i/c - Prima, in degrees; (c) Pnma - PA/nmm (in A); (e) 
PA/nmm - Pm3m (in A), where zi and Z2 refer to the two 
W-O bonds aligned predominantly with the z axis. 



ment along [001] (figure (d)). As such, all the phase 
boundaries are well-defined, although further work could 
be undertaken to locate them with greater precision. In 
the next Section we discuss the effects giving rise to this 
phase structure in greater depth. 



IV. DISCUSSION 

We see that there are two types of structural distor- 
tion occurring in WO3; rotation of the (almost rigid) 
octahedra, and displacement of the W ion from the cen- 
ter of its octahedron. The occurrence or absence of 
W offcentering is known to be determined by a bal- 
ance between electronic Coulomb repulsions (which are 
minimized for the centrosymmetric structure) and addi- 
tional bonding considerations which might stabilize the 
non-centrosymmetric phaseii. In the case of WO3, off- 
centering of the W ion results in additional hybridizations 
that lower the energy of the O 2p-like valence band com- 
pared with that of the centered structure, and raise the 
energy of the predominantly W 5d t2 g conduction band 9 . 
As the doping level is increased, and more electron den- 
sity is added to the conduction band, there is no energetic 
advantage to these additional hybridizations, and the W 
ion moves back to its centrosymmetric position. 

The fact that the small off-centering in the [110] di- 
rection is quenched first implies that the antibonding or- 
bitals corresponding to this covalent interaction are low- 
lying in the conduction band and therefore filled first. 



5 



[100] Jahn-Teller distortion 
[100] offcentring 



(a) 



(b) 



(C) 



KB B OB 
B B B B 



(d) 



(c) 



(0 



FIG. 4: Polyhedral tilts in W03, as varying with dopant 
charge. The charges in each case are, in electrons per W 
atom, are; (a) 3/16, (b) 5/16, (c) 7/16, (d) 1/2, (e) 5/8, 
(f) 1; these correspond to, respectively, the monoclinic Pc 
phase; the P2\/c phase, approaching the Pnma boundary; 
the middle of the orthorhombic phase; the Pnma - Pi/nmm 
boundary; within the stability range of the tetragonal phase; 
and the aristotype. It is clear that the distortions between 
some phases are extremely small. This poses difficulties in 
ascertaining exactly where they are located in composition 
space. 



The offcentering in the z direction persists to higher dop- 
ing concentrations and is only quenched out at ss | e~ 
per center. This 75% concentration is clearly less than 
the value of 0.98 e~ suggested by Cora et al.'s analysis 
of one-electron energies for the displacement of Re along 
[100] in ReOgS; this accords with their suggestion, how- 
ever, that crystal-field effects would favour cubic struc- 
tures, and thus any analysis under a rigid-band approxi- 
mation (or similar) would overestimate the degree of dop- 
ing necessary to cause the onset of cubicity. 

There is strong evidence of a relationship between the 
polyhedral rotation and offcentering mechanisms. In fig- 
ure0]we show sketches of the structures at different dop- 
ing concentrations to illustrate the polyhedral tilting. At 
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doping (figure^ (a)) we see that the rotations are large 
around both axes. However in figure0](b) (by which the 
offcentering in the [110] direction has been essentially 
quenched (x w yg)) we see that the rotation of poly- 
hedra in that plane also disappears, whilst rotation of 
polyhedra between connections along the z axis persists, 
as does the offcentering in that direction. This persists, 
with reducing magnitude, to x« (figure ^(c)), but by 
50% doping, (figure 01 (d)) any remaining tilt up the z 
axis (out of the page) is extremely small. However, at 
this concentration, offcentering along [001] persists. No- 
tice also that the majority of volume change in the unit 
cell, which in perovskites is caused by changing of the 
rotation angles, has also occurred by this doping level. 
Finally, any tilt systems in figures 0] (e) and (f), corre- 
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FIG. 5: Increase in energy as a function of fractional change 
in W-O bondlength for both W off-centering and Jahn-Teller 
distortion along a cubic axis. 



sponding to concentrations with centered W ions, are too 
small to observe in the figure. 

It should be noted that the very small [100]- 
offcentering - less than 0.01 A in magnitude - which per- 
sists at high doping concentration does not appear to 
be a real effect. In fact, both this and the small Jahn- 
Teller distortion at 100% doping appear to be numeri- 
cal artifacts related to the real-space mesh cutoff used: 
when the mesh is increased to 500 Ry (from 200 Ry), 
both effects disappear. As it is computationally expen- 
sive, one cannot justify this level of convergence for all 
calculations, particularly given that the order parame- 
ters delineating the phase boundaries in the system are 
well-defined at the present level of accuracy. Even so, 
it is clear that the system is very near the limits where 
Jahn-Teller distortion and/or offcentering of the W atom 
become thermodynamically stable. As expected by anal- 
ogy with ReC>3, WO3 doped with 1 electron is entirely 
cubic, with no Jahn-Teller distortion or W off-centering. 
Since the transition metal, in both cases, has a formal 
d 1 electron configuration (and d 1 ions in octahedrally co- 
ordinated complexes do not have a second-order Jahn- 
Teller distortion) the centrosymmetricity is easy to ex- 
plain. The lack of a true Jahn-Teller distortion (which 
is predicted for such a d 1 ion in octahedral coordination) 
is more subtle however, and is determined by a com- 
petition between energy lowering through hybridization 
(strongest in the cubic structure) and energy lowering 
through gap formation during the distortion. In bulk 
Re03, the broad bands and low density of states at the 
Fermi energy result in the cubic phase having the lowest 
energy. In contrast, in a molecular complex the bands are 
infinitely narrow and the Fermi energy density of states 
is infinitely large (favouring Jahn-Teller distortion). 

In figure we plot our calculated change in energy, 
as a function of the amount of both Jahn-Teller distor- 
tion and transition metal off-centering, for WO3 doped 
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FIG. 6: Mean occupancy per W atom of 5d and 6s orbitals 
versus dopant charge. 



with one electron per W center. In both cases a con- 
stant volume is maintained. The x axis shows the frac- 
tional change in bond length; for the Jahn- Teller dis- 
torted structure both W-0 bondlengths increase by this 
amount, and for the off-centered structure one increases 
and the other decreases. Note that, by this measure, 
the structure is stiffer to Jahn- Teller distortion than it 
is to off-centering. Also note that our curvature for off- 
centering is very close to that of Cora et al.'sA for Re03 
(which in turn is around half that for NaWOa). 

A Mulliken analysis of the occupancy of the three fam- 
ilies of bands due to W - 6s, 5d e , and 5d t2g - suggests 
that the occupancy of the 5d band is far from d° in un- 
doped W0 3 , and far from d 1 in the fully doped system. 
The occupancy of the d-band is closer to four in un- 
doped e-W03 (see figure |5J) and increases by approxi- 
mately 0.5 e _ in the fully doped system. The remain- 
ing half of an electron, unsurprisingly, fills the antibond- 
ing bands formed by 0(2p)-W(5d) overlap. Note that 
the mean occupancy of the 6s orbital, which is ~ 0.58 
electrons per orbital in the undoped system, hardly in- 
creases on addition of electrons. This implies a typical 
W valence population of between 4.5 and 5 electrons, in 
contrast to the idealised purely-ionic picture of a bare 
W 6+ ion (with no valence electrons); this is yet further 
evidence of a system remarkably dominated by covalent 
W-0 interactions. It is clear from further analysis that 
above 50% doping, dopant charge preferentially fills the 
5d xz and 5d yz bands; the other bands experience only 
a small change in total occupancy. This accords well 
with the suggestion that bonding-antibonding splitting 
occurs through 0(2p)-W(5c?) 7r-like overlap; the absence 
of change in xy bonding being demonstrated by the ab- 
sence in change of occupancy of the 5d xy band. However, 
given the tilting transitions and highly-deformed struc- 
tures, there is great difficulty in assigning the occupancy 
of 5d orbitals; therefore, further studies of the influence 
of band formation and orbital occupation on both first- 
and second-order Jahn- Teller distortions are ongoing. 



In terms of reconciling our calculations with previ- 
ous experimental data—, one must be somewhat cautious. 
Clarke analyzed his results under the presumption that 
there would not be significant offcentering of the W ion 
within W06 octahedra, which we believe to be incorrect: 
nevertheless, we believe that our tetragonal phase is con- 
sistent with collected data for his proposed cubic phase at 
high doping. Furthermore, our calculations deliberately 
neglect the effects of defects, chemical disorder, and other 
imperfections in the lattice. 

Finally, we mention that the original motivation of our 
study was to investigate possible self-trapping behaviour 
of charge in this material. Excess electrons in the e (low- 
temperature monoclinic) phase have been shown exper- 
imentally to self-trap, forming polarons^. We found no 
localisation of charge or deformation (polaron formation) 
in any of our calculations. This null result is of some in- 
terest, but should not be taken to mean that polarons 
do not form in this material; it is arguable that cither 
the supercell used is too small to observe polaron local- 
isation, or that the LDA will underestimate the binding 
energy of a polaron (thus causing it not to be a stable 
state of the system): further study is needed in this area. 



CONCLUSIONS 



The main conclusion that we draw from our calcula- 
tions is that the experimentally observed structural dis- 
tortions induced in Na^WOs-type bronzes by increasing 
doping are predominantly electronic in nature. By us- 
ing our methodology of adding electrons to WO3 with- 
out also adding Na atoms, we have removed any possible 
structural/disorder effect caused by the Na + cations on 
the A sites within the structure. Therefore, given that 
our calculations reproduce the experimental observed se- 
quence of symmetries upon doping ((i) Monoclinic (Pc 
and P2i/c); (ii) Orthorhombic (Pnma); (iii) Tetrago- 
nal (PA/nmm); (iv) cubic aristotype (Pm3m)) we con- 
clude that the effect of the Na + cations on the structure 
must be small. Indeed, this is not unexpected. Following 
the classification of Robin and Daji2£, as mentioned by 
Bersuker— , NaW03 is a good example of pure electronic 
doping; there is complete transfer of the donated Na elec- 
tron over to the WO3 sublattice since the 3 s band of Na 
lies well above the 5<i-antibonding band from W. 

Finally, we propose that these materials present inter- 
esting opportunities for future experimental and theo- 
retical study, given the degree of structural control that 
can be gained from doping or substitution. In particular, 
the "quenching out" of polyhedral rotation (and hence 
the opening of channels within the structure) may have 
effects on diffusion and ion intercalation in these struc- 
tures. 
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